A Prescription for Building the Milky Way's Halo 
from Disrupted Satellites 



Kathryn V. Johnston 
Institute for Advanced Study 
Olden Lane, Princeton, NJ 08540 
email: kvj@sns.ias.edu 

ABSTRACT 

We develop a semi-analytic method for determining the phase-space 
population of tidal debris along the orbit of a disrupting satellite galaxy and 
illustrate its use with a number of applications. 

We use this method to analyze Zhao's proposal that the microlensing 
events towards the Large Magellanic Cloud (LMC) might be explained by an 
appropriately placed tidal streamer, and find that his scenarios lead either to 
unacceptably high overdensities (10 - 100%) in faint star counts (apparent 
magnitudes 17.5 - 20.5) away from the Galactic plane or short timescales for 
the debris to disperse (10 8 years). 

We predict that the tidal streamers from the LMC and the Sagittarius 
dwarf galaxy currently extend over more than 2tt in azimuth along their orbits. 
Assuming that each satellite has lost half of its primordial mass, we find that 
the streamers will have overdensities in faint star counts of 10 - 100% and <1% 
respectively, and conclude that this mass loss rate is unlikely for the LMC, but 
possible for Sagittarius. If the Galaxy has accreted one hundred 10 5 — 1O 6 M 
objects (comparable to its current population of globular clusters) at distances 
of 20 - 100 kpc during its lifetime then 10% of the sky will now be covered by 
tidal streamers. 

Subject headings: Galaxy : evolution — Galaxy: Formation — Galaxy: halo - 
Galaxy: kinematics and dynamics 
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1. Introduction 

Two fundamental astrophysical questions that can be addressed with Galactic 
research are: how are galaxies formed?; and what are they made from? If galaxies formed 
hierarchically, we expect satellite accretion to have played some role in Galactic history, and 
might hope to address the former question by searching for signatures of these events. The 
latter question is currently being addressed by microlensing surveys of the Galactic halo 
(e.g. Alcock et al. 1996|) . In most models for the currently available microlensing results 



it is assumed that the halo is sufficiently smooth for its global baryonic population to be 
determined from a local sample. |Metcalf fc Silk (1996)| and |Zhao (1997)| have pointed out 



that lumps in the halo, either caused by dark clusters or debris from satellite disruption, 
would jeopardize this interpretation. 

To address both these questions, therefore, we need to quantify the significance of 
the presence of accreted structure in the Galactic halo. [Unavane, Wyse fc Gilmore (1996 



looked at the distribution of halo stars in age and metallicity and concluded that no more 
than 10% of these stars could have been been accreted from the destruction of galaxies 
with stellar populations like those of the dwarf spheroidal satellites of the Milky Way. An 
alternative approach, which requires no assumptions about the age or metallicity of stars 
in accreted satellites, is to look for signatures of accretion in the kinematic and spatial 
distribution of the halo. 

Numerical simulations have confirmed the physical intuition that lumps observed in 
the halo's phase-space distribution could be associated with accretion events. They have 
shown that an initially spherical satellite will become distorted when introduced into a 
tidal field, with stripped stars forming long streams along its orbit, both ahead and behind 
(e.g. McGlynn 19901 ; |Moore fc Davis 1994| ; |(Jh, Lin fc Aarseth 1995] ; |Piatek fc Pryor 1995 



|Velazquez fc White 1995| ; |Johnston, Spergel fc Hernquist 1995] ; Johnston, Hernquist & Bolte 



1996 — hereafter XFTB; Kroupa 1997| ). Analyses of the tidal debris in these simulations have 



demonstrated that these star streams can maintain coherence for the lifetime of the Galaxy, 
and that long-lived moving groups of stars are a natural consequence of tidal disruption 
QJohnston, Spergel fc Hernquist 1995] ; |JHB| ) . 

There is increasing evidence of morphological distortions, reminiscent of these 
simulations, in the Milky Way's satellites: surface density contours of the Sagittarius dwarf 
galaxy (discovered by |lbata, Gilmore fc Irwin 1994| ) have now been reported to cover an 
angular region 8 x 22 square-degrees QMateo et al. 19961 ; |Alard 19~96l ; [Fahlman et al. 1996| ; 



Ibata et al. 1996j ); distorted isophotes have been seen in several globular clusters (|Grillmair| 



et al. 1995]) ; and "extra-tidal" stars have been discovered in most of the dwarf spheroidal 



satellites ( Gould et al. 19921 ; |Irwin fc Hatzidimitriou 1995| ; [Kuhn, Smith fc Hawley 1996Q . 
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Substructure in the spatial distribution of the Milky Way's dwarf spheroidal satellites, 
in the form of alignments along Great Circles, has also historically been attributed to the 
ancient disruption of a much larger body ( [Lynden-Bell 1976] ; |Kunkel fc Demers (1977)D , 
and further possible associations among the halo's globular cluster population have recently 
been found (|Lin fc Richer 1992| ; Majewski 1991 ; Fusi-Pecci et al. 19951 ; [Lynden-Beir&] 



Lynden-Bell 1995|) 



In the stellar halo, the advent of various digitized Plate Surveys (e.g. Maddox~et| 
al. 1990| ; Pennington et al. 1993| ; |Weir, Fayyad fc Djorgovski 1995|) raises the possibility 



of searching for lumps in its projected angular distribution: the dwarf spheroidal galaxy 
Sextans was found as a simple overdensity in star counts in the APM survey ( |lrwin et al.| 
1990| ); [Klenya et al. (1997)| have developed a circular filter technique that can be used to 



search for systems with similar morphologies and produce a quantitative estimate of the 
completeness of the sample that has been found; and motivated by the work on satellite 
alignments |JHB| proposed the method of Great Circle Cell Counts to look for streams of 
stars along Great Circles that might result from the disruption of a satellite galaxy. 

Observers have been aware of moving groups as distinct lumps in the local stellar 
phase-space distribution for a long time (e.g. |Eggen 19U5| ) and similar structures have also 



now been reported in the halo (e.g. [Sommer-Larsen fc Christiansen 1987| ; |Arnold fc Gilmore] 



1992| ; [Majewski, Munn fc Hawley 1994] , 1996; — for a comprehensive review of this subject 



sec 



Majewski, Hawley fc Munn 1996|) . These discoveries are currently being followed up 



with systematic spatial and kinematic surveys of larger portions of the halo with the hope 
of assessing how exceptional this substructure is (Majewski - private communication). 

Overall, the qualitative comparison of numerical and observational work outlined above 
clearly demonstrates that satellite accretion is ongoing in the Milky Way, and is likely to 
have occurred often in the past. To address the aim of quantifying the presence of accreted 
structure, the characteristics of a specific phase-space feature (such as a moving group) or 
the frequency with which such features are found in the stellar halo could be interpreted 
using a suite of N-body simulations. However, the number of accretion events that such a 
study can survey is limited by time and computational cost considerations. An alternative 
approach was adopted by Tremaine (1993), who restricted the description of a trail to 
simple estimates for its length and width given the parent satellite's mass and distance 
from the Galactic center. This allowed him to derive analytic expressions for the fractional 
sky-coverage of tidal debris from an ensemble of satellites, but at the expense of a more 
precise knowledge of the structure of the individual trails. 



In this paper, we develop a semi-analytic technique for finding the density distribution 
along the entire length of a debris trail, given a satellite of any mass and orbit, assuming 
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only a disruption timescale. This method is complimentary to both previous numerical and 
analytic work, being less computational expensive than the former and hence able to cover 
a wider region of parameter space, while providing a more detailed description of individual 
trails than the latter. The success of the method is tested through comparison with N-body 
simulations. As an example of its use we discuss the implications of this model for the 
surface density structure of debris trails in the Milky Way system. In a further paper we 
will extend the technique to generate full phase-space models of halos with a variety of 
accretion histories, and examine their observable properties flJohnston fc Majewski 19971) . 

We present the method in §2, apply it to the Milky Way system in §3, and use it to test 
the plausibility of Zhao's (1997) model of the microlensing results towards the LMC in §4. 
We summarize and discuss the limitations of this approach in §5. The N-body simulations 
used to confirm the validity of our technique are outlined in Appendix A. 

2. Method 

In this section we present our method for constructing the phase-space structure along a 
debris trail. We outline our understanding of the debris' orbital energy distribution, gained 
from studies of tidally disrupting hydrodynamic systems, in §2.1. In §2.2 we apply this 
picture to find the energy distribution of tidal debris from stellar systems, and then convert 
this to a density distribution along the tidal streamers by using the general properties of 
orbits in a halo potential. This approach is tested by comparing its predictions with the 
final particle positions in N-body simulations which have followed satellites evolving in a 
Milky Way potential for 10 Gyrs (see Appendix A). We summarize the technique in §2.3. 

2.1. Hydrodynamic Systems 

Analytic descriptions of tidal disruption and debris dispersal have been developed for 
hydrodynamic systems such as planetesimals encountering planets in the Solar System (e.g. 
|Sridhar fc Tremaine 1992] ) and stars falling into black holes (e.g. [Evans fc Kochanek 1989| ; 
Kochanek 1994j ). In these studies, the satellite is assumed to move on a parabolic orbit, 
with internal properties such that it disrupts at pericenter. The size of the satellite as it 
disrupts is the physical scale where the internal and tidal forces balance — roughly 



/ m \ 1/3 

rtide = UJ 



(1) 



- 5 - 



where m is the mass of the satellite system, R is the distance of closest approach between 
the parent and satellite and M is the mass of the parent system (cf. |Roche 1847| ). The 
typical scale of the energy distribution of the debris is then the change in the Keplerian 
potential energy on this spatial scale, 

# _ /m\l GM 

eKepler = T tide — = ^ — j — . (2) 

This energy scale lies between the orbital energy, e or b = GM/R = (M/m) 1//3 eKe P ier and 
the internal energy of the satellite, e- mt = Gm/r = (m/M) 2 ^ 3 eK cp icr- Thus, in the regime 
typically considered, where (m/M) 1 ' 3 <C 1, the debris from a disruption event occupies a 
region in energy space where the perturbation from the satellite's original orbital energy is 
small compared to e or b, but large compared to e- m t. 

Indeed, Sridhar & Tremaine (1992) and Evans Sz Kochanek (1989) both found the 
energy distribution of the debris to be continuous over a finite range, with scale eKopior, 
centered on the initial orbital energy. Approximately half of the material moved to 
unbound, hyperbolic orbits and half moved to bound elliptic orbits. Ignoring the influences 
of self-gravity and fluid effects, the subsequent evolution of the debris could be followed 
analytically by simply mapping the distribution in energy to one in orbital time-periods. 



2.2. Stellar Systems 

The same physical principles used to determine the energy distribution and dispersal 
of tidal debris for hydrodynamic systems apply to stellar systems, but the theory is less 
complex since it involves only collisionless particle dynamics. Its application to Galactic 
satellites differs in two ways: (i) the parent potential is not Keplerian; and (ii) the satellite 
undergoes repeated encounters. Mass is lost continuously, but predominantly during the 



pericentric passages (see e.g. |JHB| , their figure 2). The implications of these differences are 
outlined in §2.2.1 and §2.2.2, and the resulting prediction for the matter distribution along 
a tidal streamer is presented in §2.2.3. The effect of the omission of self-gravity from this 
analysis is discussed in §2.2.4. 



2.2.1. Orbits in a Halo Potential 

Orbits in a spherical potential can be classified by their energy E and angular 
momentum J, which completely specify their planar paths, (R, \l/) in polar coordinates in 
the orbital plane, as well as their radial and azimuthal time periods (Tr,T^). Hence, for 
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the purpose of characterizing the properties of orbits in the halo we consider a spherical 
approximation to the full Galactic potential that is used in the numerical simulations (see 
Appendix A). This potential $mw is constructed from bulge and halo components taken 



directly from equation (Al|) and a third spherical component $ = — GM^^/ \/R? + a 2 , with 



the mass and scale length of the disk given in equation (|A1|) . 




For this spherical potential, the extent of the energy distribution in the debris from a 
satellite on orbit with pericenter R p can be found analytically by using generalized forms of 
equations ([[J) and @ 

' m ' ~ (3) 

£ * = r -35=U) ^5S = (Gm ^» / (dSj ' (4) 

where M p = Rp(dQ/dR)/G is the mass of the Galaxy within R p []. However, the azimuthal 
and radial orbital time periods must be found numerically. 

The upper panel of Figure ^ show T# as a function of E for J / J circ = 0.05, 0.1, 0.15, 1.0 
(where J C irc is the angular momentum of a circular orbit). Note that the lines in this panel 
are almost indistinguishable as depends only very weakly on J. For a circular orbit 
can be calculated exactly, 

The solid lines in the lower panel of Figure p] show the ratio Tr/T^ for the same values of 
J I Jcirc- For a near-circular orbit, Tr was calculated using the epicyclic approximation (see 
[Binney fc Tremaine 1987| ) 



Tdrc _27T K 2_i± ,SJ^_d^_ 3*& 
R ~ «' dR* + R* ~ dR? RdR W 

For a radial orbit Tr/T^ = 1/2. 

Since we are primarily interested in the construction of the outer halo (> 20 kpc) from 
tidal debris, we might hope to represent the full Milky Way potential as logarithmic in this 
region, 

$i og = v^ ivc log(R/Ro), v circ = 200 km/s, Ro = arbitrary constant. (7) 

(The I.A.U standard for the circular velocity at the solar radius is Vq =220 km/s, but recent 
estimates tend to be lower than this — see [Sackett 1997| for a summary.) This further 



1 In a calculation analogous to Roche's (1847) work, King (1962) found the limiting radius of a star cluster 



on an orbit of eccentricity e and pericenter R p to be rtidc = R P (m/ M (3 + e)) 1 / 3 . 
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approximation is not strictly necessary, but it allows us to express equations (||), @ and 
([]) in simple form: 






(9) 



"'circ 



R 



'circ 




(10) 



In Figure |^ we make some assessment of how successful this representation is for our 
own three-component model of the Milky Way. The dotted line shows (ei og /eMw — 1) as a 
function of R p , and the solid line shows (T log /T M w ~ 1) f° r time periods of circular orbits 
as a function of -R C irc, where "log" and "MW" denote quantities calculated in the potentials 
$i og and $mw respectively. The Figure demonstrates that these properties in the two 
potentials differ by only 10% in the region of interest. 

We use Figures |l| and || to justify the following approximations: 
1. We use the logarithmic potential to estimate the scale of the energy distribution (eq. 



2. We assume that the azimuthal time period Tq, of an orbit is a function only of its 
energy, independent of its angular momentum or eccentricity. We calculate using 
the expression for a circular orbit of energy E in a logarithmic potential given by 
equation @. This allows us to write down the relation between T# for orbits of 
energy E and (E + AE) as 



3. We take Tr/T^ = constant = yl/2. Assuming this condition is strictly true, then the 



angular difference between successive turning points (or pericenters and apocenters) 
is the same along all orbits. Given the azimuthal time dependence ^(t) of an orbit 
of energy E and period T^, we can express the azimuthal path of one with energy 
(E + AE) and period tT# as (\l/ + ^(*o + V r ))> where \l/ an d t are constants chosen 
to take into account the difference between the azimuthal and radial phases of the 
two orbits at t — 0. 




(ii) 
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2.2.2. Energy Distribution in Debris 

To simplify the calculation of the energy distribution in the tidal streamers we model 
the disruption of a satellite as a series of discrete mass loss events occurring on each 
pericentric passage (which will clearly be most appropriate for eccentric orbits). 

Figure |3] is a "snapshot" from a simulated encounter, showing the positions of particles 
in the orbital plane of a satellite. The dots in the plot show particles still bound to the 
satellite and the bold circle shows the extent r tide from equation ([$]) that sets the scale of 
the energy changes. The dashed line shows the satellite's orbit and the plus (+) and minus 
(— ) signs show unbound particles which have moved to orbits of higher and lower energies 
respectively. The plot suggests that mass lost on passages prior to disruption will make 
two distinct contributions to the energy distribution in the tidal streamers, at positive and 
negative AE, clearly separated from each other in energy by the region where AE is too 
small for a particle to escape from the satellite. In contrast, in their study of hydrodynamic 
systems disrupting during a single encounter, Evans & Kochanek (1989) found the density 
of debris to be constant in energy across the energy scale eKepier centered on the satellite's 
orbit (i.e. dN/d(AE) = l/2e for — e < AE < e). In the light of this study, we expect the 
particles which become unbound on the final passage, as the satellite disrupts completely, 
to remain at energies close to the original orbit rather than forming distinct populations at 
positive and negative AE. 

To convert these expectations to a formalism for the finding the full energy distribution 
of mass lost from a satellite that disrupts over several orbits, we first define the energy scale 
of debris lost on the ith of n pericentric passages to be £j = (Gmj/i? p ) 1 / 3 (^ irc ) 2//3 , where 
rrii is the mass still bound to the satellite on that passage (cf. eq ||). Let iV peri (g i ) and 
-^disrupt (<?«) be the fractional number densities in scaled energy = AE/ei for particles 
unbound on the ith pericentric passage and as the satellite disrupts respectively. The full 
distribution of tidal debris is given by 



d{AE) ~ J- 



Srrii N pcri (AE/ej 

^■debris 



| m n N disTupt (AE/e n ) 

^■debris 



where mdebris is the total mass of the debris and Srrii = (rrii — m i+ i) is the mass lost on the 
ith passage. 

Figure |] shows the number density dN/dq of unbound particles in scaled energy 
changes q = AE/e at the end of each of the simulations (histograms - see Appendix A) 
where e is calculated from equation fl8|) with the values for m and R p listed in Table [I]. 
As predicted, the debris in all the models covers a similar range in q despite the variety 
of satellite orbits and masses employed. In the models that were not disrupted during 
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Model 


m 


R p 


R a j Rp 


n 


/ 




(M ) 


(kpc) 








1 


1.0 x 10 9 


39.6 


3.4 


5 


0.18 


2 


4.1 x 10 7 


43.4 


4.0 


4 


0.30 


3 


1.9 x 10 7 


45.2 


2.4 


1 


0.65 


4 


6.9 x 10 7 


26.4 


3.0 


5 


0.30 


5 


3.3 x 10 s 


43.4 


4.0 


4 


0.30 


6 


1.2 x 10 8 


43.4 


4.0 


4 


0.30 


7 


1.4 x 10 7 


43.4 


4.0 


4 


0.30 


8 


5.1 x 10 6 


43.4 


4.0 


4 


0.30 


9 


5.8 x 10 7 


43.4 


4.0 


2 


0.50 


10 


3.3 x 10 7 


43.4 


4.0 


4 


0.25 


11 


2.9 x 10 7 


43.4 


4.0 


4 


0.20 


12 


2.4 x 10 7 


43.4 


4.0 


4 


0.12 



Table 1: Parameters of models used to generate debris trails. Columns: (1) model number; 
(2) satellite mass; (3) pericenter of orbit; (4) ratio of apocentric to pericentric distance; (5) 
number of pericentric passages either prior to destruction, or during simulation; (6) fraction 
of bound mass lost on each passage. 
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the course of the simulation (Models 1, 10, 11 and 12), the debris shows two distinct 
populations in energy changes, as expected from the discussion above. The other models 
show similar distributions at lower and higher energies formed at pericentric passages prior 
to the destruction of the satellite, along with a central population of particles with smaller 
\AE\ from the final disruption event. The histograms are similar enough to suggest that 
we can postulate forms for intrinsic distributions in q for the material torn during each 
pericentric passage as 



1 



(\q\ — a)/(b — a), if a < \q\ < b; 



N P e ^ = — a X \ ( c ~ l<?D/( c - 6 )> ifb<\q\<c; (13) 

0, otherwise, 

and from the satellite disruption as 

W!)44l' if i g ' <rf; ( 14 ) 
2d [ 0, otherwise, 

where a = 0.6, b = 1.1, c = 3.1 and d = 1.5. The shape of these distributions and values for 
the parameters are chosen by simple inspection of the histograms in Figure 

To verify the proposed forms for the scaled energy distributions with each simulation 
we need to specify the number of pericentric passages n and the mass of the satellite on each 
passage m*, and use these to convert iVp OT i and Abrupt given by equations (0) and (|14D 
to the full distribution in energy changes dN/d(AE) using equation (|T2D . The mass loss 
history — which we characterize by the values of — depends non-trivially on the mass, 
physical scale, and orbit of the satellite, and should be investigated with simulations if a 
detailed model of a specific tidal disruption event is needed. For our purposes of generating 
physically reasonable debris trails, we assume that each satellite loses a constant fraction 
/ of its remaining mass instantaneously on each passage. To compare with the simulations 
we choose / and n to approximate their mass loss histories by visual inspection of figure 2 
in |JHB| (see Table p. Substituting Srrii = frrii and m.j = (1 — f) l m in equation ([12]) we find 
the dN/d(AE) plotted as bold lines in Figure £|. The fits are good enough to justify our use 
of the simple functional forms for the energy distributions given by equations [13] and [TJ|. 



2.2.3. Matter Distribution Along Tidal Streamers 

For any given orbit, with azimuthal time dependence *&(t), we can generate a detailed 
description of the density structure along the tidal streamers at any time by combining 
the energy distributions and mass loss histories outlined in §2.2.2 with the approximations 
for orbital properties given in §2.2.1. We assume that particles lost at pericenter initially 
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all have the same radial phase as the satellite. Then the problem of finding the number 
density of particles at ^ from mass lost with energy distribution dN/d(AE) during the ith 
pericentric passage at time tj reduces to a mapping (using eqs. [§] , |TT| and 

dNi dNi d(AE) dr 

~d¥ ~ d(AE) dr d^(U + (t - U) /r) 
v 2 1 



v 2 

circ 




Gmi/R p J ' u^(U + u) 



(15) 



where u = {t — ti)/r, qi = AEjti and N(qi) is given by either equation (|i~3|) or (|14l), as 
appropriate. The complete distribution is calculated by summing over mass lost on each of 
the pericentric passages. 

Figure || shows the fraction of unbound particles per degree, dN/d^! , along the orbit of 
the satellite at the end of each of the simulations (histograms). The bold lines in the plot 
are the predictions of the mapping given in equation (|15D, using the parameters listed in 
Table [I] and taking \1/ to be the angular velocity along the (non-planar) path of the satellite 
found by integrating its orbit in the full Milky Way potential given in Appendix A. In 
general, the semi-analytic approach reproduces the density distribution to within a factor 
of two and over two orders of magnitude in amplitude. 



2.2.4- Self- Gravity in Debris 

The gravitational interactions of particles in the debris trails are neglected in both the 
simulations (since the code uses spherical harmonic expansions to evaluate the potential and 
cannot represent the geometry of the streamers) and the above semi-analytic description. 
However, we can argue that self-gravity is less important than tidal forces for the predicted 
angular distribution by comparing the influence of the two at each point along the streamers. 



Following |Roche 1847] and |King (1962)| , we expect the internal and external forces to 



balance when the debris has a density of order 



Pcrit = — ft3~- ( 16 ) 



The debris density p can be calculated from equation ( |15D by including estimates for the 
stream's angular width w and radial extent AR. The former is fixed by the range in orbital 
inclinations, set at pericenter to 

w = (r tidc /R p ), (17) 
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and the latter is determined by the scale of the energy changes, A_R = e/(d$/dR) = Rw 
(cf. eq. 0). Hence 



p(R, ) 



m debris dN m d ebri S f R p \ 2 dN 



R 2 wAR d^f R 3 \r tid J d^ y 

For the center of a disrupting system, the ratio of this density to the critical density for 
self-gravity to be important is 

/? I i 

(19) 



p _ R p 1 1 



p cr it R 2dm 
(combining eqs. [TJ, |5|, @) and @). 



This ratio is plotted in Figure^ for orbits with J / J circ = 0.28, 0.36... 1.0 (or eccentricities 
R p / R a ~ 0.07.. .1.0, where R a is the apocentric distance). At pericenter the curves reach 
minima whose amplitudes decrease with J as ^(R p ) increases. At apocenter there is a 
corresponding decrease in ^(R a ) with J, but this is offset by an increase in the cross 
sectional area of the trail because at large distances it is less confined by the tidal field. 
Thus the density evolution along a circular orbit (shown in bold) serves as an upper limit 
to the behavior along any orbit. Since this never exceeds the critical density we conclude 
that our prediction for the matter distribution along the streamers is consistent with the 
neglect of self-gravity in the analysis. (The density becomes singular as t —>■ since the 
initial angular extent of the debris in our description is effectively zero at this time.) 

Note that mass escaping over the tidal radius of a still-bound satellite will also have 



density <; p crit , and is spread over a similar energy scale (cf. eqns jnj and JHJ). Hence our 
neglect of self-gravity is equally valid when modeling the behavior of tidal streamers from 
Milky Way satellites (such as the LMC). 



2.3. Summary of Technique 

Our method for generating debris distributions along tidal streamers, can be 
summarized as follows: 

(i) We assume that self-gravity is negligible in the tidal streamers. 

(ii) We apply a simple physical picture to understand the energy scale, e, around the 
satellite's orbital energy over which the debris is distributed (see §2.1, §2.2.1 and 
eq. p[) and take a functional form for the full distribution in scaled energy changes, 
N(AE/e), by examining the results from N-body simulations (see §2.2.2 and eqs. |H| 
and 0) 
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(iii) We make the approximation that orbital time periods depend only on their energies. 
Using this approximation, along with the characteristics of orbits in a purely 
logarithmic potential (eq. J7|), we relate the energy change AE to a change in orbital 



time period (§2.2.1 and eq. |ll|). 



(iv) We model mass loss as occurring instantaneously at the satellite's pericentric passages 
and assume that the subsequent turning points along the debris' and satellite's orbit 
are aligned in azimuth (or Tr/T^, =constant). 

(v) Combining the above points with the satellite's orbit (R(t), ^(t)), we predict the debris 

distribution along the path of the satellite from mass lost at each previous pericentric 
passage, and sum these to find the full density profile (§2.2.3 and §2.2.4). 



3. Applications to the Milky Way System 

We illustrate some uses of the technique developed in the previous section by asking: 
is it possible to detect tidal streamers from any of the Galactic satellites? (§3.1); how long 
will overdensities in star counts be noticeable following a disruption event? (§3.2); what 
is the angular scale of a tidal streamer and what fraction of the sky might an ensemble of 
such streamers cover? (§3.3). 



3.1. Should We Expect to See Tidal Streamers from Galactic Satellites? 

If we know the position, radial velocity and proper motion of a satellite, we can 
integrate its orbit backwards in time in our (or any other) Milky Way potential to find 
its angular path ^(t), pericentric position R p and number of encounters n in the last 10 
Gyrs. Assuming the satellite, currently mass m sat , has lost mass mdebris in this time (or, 
equivalently, a fraction / = [mdebris/ (m de bris + m. sat )] 1//n of its bound mass on each passage) 
we can use equation flT5|) to calculate the matter distribution dN/ along the path of the 
satellite. This can be normalized to give a mass surface density in M /(degree) 2 

v-i/.t.\ m dehris dN R p dN 

^ ) = —M =md ^^- c M (20) 



where w is the width of the debris trail (see eq. ||17| ). (For simplicity, we have adopted the 



viewpoint of an observer at the Galactic center in this and all subsequent calculations.) 

Figure |7| shows the prediction for £ along the paths (calculated from their proper- 
motions, see references in Table 0) of the LMC and the Sagittarius, Ursa Minor and Draco 
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Satellite 


mass 


L 


R 


Ra/Rp 




debris scale 


proper 




(M Q ) 


(L Q ) 


(kpc) 




(Gyrs) 


(width x length) 


motion? 


LMC 


10 10 


2 x 10 9 


31.9 


2.2 


1.5 


17° x 360° 


a 


Sgr 


10 8 


2 x 10 7 


12.7 


3.1 


0.7 


6° x 360° 


b 


Scul 


1.4 x 10 7 


1.6 x 10 6 


52.0 


2.4 


2.6 


2° x 36° 


c 


UMinor 


3.9 x 10 7 


3.0 x 10 5 


45.7 


2.0 


2.0 


3° x 90° 


d 


Draco 


5.2 x 10 7 


2.5 x 10 5 


76.0 




2.6 


3° x 60° 




Carina 


1.1 x 10 7 


2.9 x 10 5 


89.0 




3.0 


2° x 30° 




Sextans 


2.6 x 10 7 


8.3 x 10 5 


91.0 




3.1 


2° x 45° 




Fornax 


1.2 x 10 8 


2.5 x 10 7 


133.0 




4.5 


3° x 36° 





Table 2: Parameters of Milky Way satellites. Columns: (1) name; (2) mass; (3) luminosity; 
(4) distance from Galactic center; (5) ratio of apocentric to pericentric distance; (6) azimuthal 
time period; (7) approximate debris scale, calculated from the width and length estimates 
given by equations fll7|) and ( p5|) ; (8) proper motion references: a. IJones, Klemola fc Lin 
19941 ; b. |lrwin et al. (1996)1 ; c. [Schweitzer, Cudworth, Ma.jewski fc Suntzeff (1995)1 ; d. 



Schweitzer, Cudworth fc Majewski (1997)| . Columns (4) and (5) give pericentric distances and 
time periods integrated in the Galactic potential for those with proper motions. Otherwise 
these give the current distance of the satellite and the time period for a circular orbit at 
that distance. The mass, luminosity and current galactocentric distances are from Jones, 
Klemola & Lin (1994 — LMC), Ibata, Irwin & Gilmore (1994 — Sgr) and Mateo et al. (1993 
— remaining dwarf spheroidals) . 
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dwarf spheroidals assuming mdebris = ms&t f° r each satellite. From estimates for their 
central luminosity densities and mass-to-light ratios, (see, e.g. [Irwin fc Hatzidimitriou 



the central surface densities of these satellites are typically > 10 8 M Q /degree 2 . The 
current position of each satellite is given in Galactic latitude and longitude (l,b), and the 
corresponding value for the position \l/ along the orbit is indicated with the vertical dotted 
lines. Note that since the extent of the debris' energy distribution is only weakly dependent 

1/3 

on the instantaneous mass of the satellite (e oc ) the shape of these curves is nearly 
independent of mdebris so long as mdebris < m sat , and they can be simply scaled down in 
amplitude for lower mass loss rates. 

Debris from the LMC and Sagittarius spreads over more than 2n in \1/ in 10 Gyrs, and 
in these cases the plot shows the result of summing over all mass which has the same angular 
phase (i.e. assuming that the orbit is perfectly planar and the same phase corresponds to 
the same line of sight). Our prediction for the extent of the debris trail from the LMC is 
consistent with the hypothesis that an ancient accretion event is responsible for the curious 
alignment of several dwarf spheroidal satellites ( Lynden-Bell 1976| , 1982), young globular 
clusters ( [Lin fc Richer 1992 ; Majcwski 1994 ; Fusi-Pecci et al. 1995 ) and proper motions of 



some of these objects ( [Majewski, Phelps fc Rich 1996|) along the Great Circle define by the 



LMC and the Magellanic Stream. This calculation also confirms the idea that a globular 
cluster that was previously a member of Sagittarius (and several are currently associated 
with this dwarf — see pa Costa fc Armandroff 1995| ) could now be on the opposite side of 
the Galaxy ( [Fusi-Pecci et al. 1995|) since the debris is expected to encircle the sky within 



the lifetime of the Milky Way. 

For any given luminosity function, we can use our estimate for E to find the surface 
density of stars along the trail in apparent magnitude bins. We illustrate this idea by 
considering horizontal branch (HB) stars, which we expect to be useful tracers of debris 
trails. Although not as bright as giant stars, they are easier to distinguish from the 
foreground disk population using their blue colors. Moreover, we can approximate their 
absolute magnitudes as 

M HB = 0.66 + 0.19[Fe/H] (21) 



( |Walkcr 1992| ), and therefore find their distances. From equation (ET]), we know that the 



apparent magnitudes of HB stars at distances of 25 - 100 kpc lie in the range 17.5 - 20.5. In 
this example we will ignore any color information, and simply contrast their number density 
along each streamer with star counts at the same Galactic latitude and longitude (I, b) and 
in the same apparent magnitude range, estimated from the Bahcall-Soneira Galactic model, 
iVes flBahcall fc Soneira 1980j ). We convert Figure [5] into an estimate for the number of 
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satellite horizontal-branch stars per square degree via 

iV HB = n HB — (22) 

m sat 

where tt-hb is the number of horizontal branch stars per unit luminosity for the satellite's 
population, and L sat is the total luminosity of the satellite. To estimate a typical value 
for nnB we note that the observed ratio of the number of HB stars to RGB stars with 
magnitudes less than the horizontal branch lies in the range n^B /^rgb ~ 1 — 2 ( |Mighell 
1990| ; |Buonanno Corsi fc Fusi Pecci 1985|) . Then, using [Bergbusch fc VandenBerg's (1992) 



isochrones for an 8.0 Gyr, [Fe/H] = —1.03 cluster as a representative population, we find 
n RGB = 1/528, which implies n HB ~ 0.002 - 0.004. 

In Figure || we plot A^b/A^bs f° r ^hb = 0.003. The squares mark positions along each 
satellite's orbit that are less than 20 degrees from the Galactic plane, where the background 
counts are largest. The comparison suggests that even for large mass loss rates, only the 
tidal streamers from the LMC would be detectable as local over-densities in star counts, 
though the debris from the other satellites might be found by integrating the star counts 
along the suspected orbital path (as proposed in JHB) or using additional color, kinematic 
or distance information. These examples also imply that tidal streamers from such minor 
accretion events in external galaxies are likely to be unobservable. 



3.2. How Long Will We See a Satellite After Disruption? 

The central surface density of a remnant at time t after the destruction of a satellite is 
given by 

So sat cN/r tidc £ 2d [ ] 

where S = m sat Rp/r^ ide is the approximate surface density of the satellite as it disrupts. 
This is plotted in Figure [] for a circular orbit (bold) and orbits with eccentricities spanning 
the range R v j R a = 0.07 — 1.0, where R a is the apocentric distance (as in Figure ||). In all 
cases, E/E oscillates between apocenter and pericenter, but with a characteristic drop in 
amplitude by more than a factor of 10 within the first two radial periods (~ \/2T$) — in 
agreement with the results from Kroupa's (1997) simulations. For objects that have already 
been destroyed, we may take the surface densities of the dwarf spheroidals observed today 
(which are 10-100 times the background star counts — see [Irwin fc Hatzidimitriou 1995| ) 



as upper limit to S . This implies that the signatures of such events as a local overdensity 
in star counts would be clear only for the time period of the orbit (i.e. ~ 1 Gyr, see eq. 
0). However, the amplitude of the oscillation only decays linearly with time. This suggests 
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that either a matched filtering technique, where star counts are integrated over some region 
(e.g. |J Hrj| ) , or a search including kinematic information could effectively find debris for the 
lifetime of the Galaxy (tens of orbital periods). 



3.3. How Lumpy is the Halo? 

We can estimate a timescale to characterize the rate of debris dispersal by calculating 
how long it will take the streamers from a satellite to extend 27r along its orbit: 

Zk = = \ = 1 m 

X* r(2e)-l l-exp(-e/v£) 1 - exp(-(2Gw/iV4c) 1/3 ) ' ^ > 

Figure |10| plots T 27r /T^ as a function of pericentric radius for various satellite masses, and 
for each of the Milky Way's satellites (at either their pericentric or current positions for 
those with or without proper motion data respectively). The upper axis shows T# for a 
circular orbit at radius R p . The plot confirms the result of §3.1 that debris from the LMC 
and Sagittarius is expected to spread more than 2tt along their orbits over the lifetime of 
the Galaxy. 

We can use this timescale to find the approximate angular length L of a streamer from 
a satellite that has lost mass during the last 10 Gyrs: 

L = ^=V, (25) 

and this is included in Table 0. Combining this with the width of the streamer (eq. [ |17|j ) 
and approximating T# with the time period of a circular orbit at the pericentric position, 
we can estimate the sky coverage w x L of debris from satellites of various masses and 



pericentric distances, as shown in Figure |TT 



The lumpiness of halos with different accretion histories can be assessed by 
superimposing an ensemble of these trails and calculating how much of the sky they cover. 
We assume that M acc of the mass of the halo has been accreted in the form of M acc /m sat 
satellites of mass m sat , to form a density distribution that follows that of the stellar halo 
(where p oc r -3 ' 5 , see [Freeman 1996| ). Then we integrate the sky coverage of each individual 
trail over the number density of accreted satellites to find the expected number of trails 
^trails intersecting any random line of sight (or, equivalently, how much of the sky is 
covered) . 



Figure [L2| shows ra tra ii s as a function of satellite mass for M acc = 10 , 10 , 10 , 10 M & 



see also [lremaine 1993| for a simple — and elegant — version of this calculation). Of 
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course, we know neither the accretion history of the Galaxy nor the primordial spectrum 
of satellite masses. However, semi-analytic calculations suggest that many of the Galactic 
globular clusters will disrupt in the next Hubble time and that those we see today are the 
remnants of a much larger primordial population (e.g. Murali fc Weinberg 1996j ; Gncdin Ec 



striker 19971 ). From §3.2, we would not expect to see trails from such disruption events 
for more than a few orbital periods in star counts alone. However, if we suppose that 
~ 100 similar objects (with masses of 10 5 - 10 6 M Q ) have been accreted in the halo in the 
last Hubble time, then Figure [12| predicts that their debris trails would cover 10% of the 
sky. This significant sky coverage encourages the idea that many could be detected with 
kinematic surveys of halo stars along random lines of sight. 



4. Microlensing in a Lumpy Halo 



In a recent paper, [Zhao (1997)| proposed several different scenarios in which the 



superposition of a debris trail between us and the Magellanic Clouds could explain all 
the microlensing events observed by the MACHO team ( [Alcock et al. 1996| ). In a field 
2° northwest of the LMC bar, [Zaritsky fc Lin (1997)| find tantalizing evidence for such a 



foreground population of a suitable density and distance: 5% of the clump giants in their 
survery are 0.9 magnitudes brighter than the mean (note that Gould 1997| questions the 



conclusion that this population could indeed account for the microlensing events). With the 
tools developed in this paper we can re-examine the likelihood of each of Zhao's scenarios, 
and hence the interpretation of Zaritsky Sz Lin's clump stars as tidal debris. 

In Zhao's first scenario, he assumes that the Magellanic plane (defined by the LMC 
and the Magellanic stream) is populated by a stream of debris, with lensing mass 1O 1O M , 
uniformly distributed along a great circle — equivalent to the model we discussed in §3.1 
for tidal streamers from the LMC. Our formalism predicts that the streamers would cover 
an area 17 x 360 square-degrees (see Table |2|). If the debris has the same mass-to-light ratio 
as the LMC, this would correspond to the luminosity of the entire stellar halo (~ 1O 9 L , 



see [Freeman 1996|) in only 10% of the sky and would produce an over-density in star-counts 
across the Magellanic Plane at high Galactic latitude of 10-100% (see Figure §). It seems 
unlikely that such a significant population could have been missed in previous studies since 
the number-counts predicted by the Bahcall-Soneira model have been found to provide a 
reasonable fit to observations in all fields to which they have been applied ( [Bahcall 1986j ; 



Reid fc Majewski l993j). 



As an alternative, Zhao suggests debris just covering the extent of the clouds, either 
1O 1O M immediately in front of the clouds (~ 45 kpc from the Galactic center), or 
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10 8 — 1O 9 M at the Galactocentric distance of Sagittarius (16 kpc). In both cases, the 
timescale for such debris to disperse over an orbit is ~ 3Gyrs (from the positions of the 
LMC and Sagittarius in Figure [II]) which suggests that such an alignment over just 10° 
would only last for < 10 8 years. Thus, this requires the destruction of an object of either 
mass within the last 10 8 years, and the alignment of its debris in just l/400th of the sky. 
We conclude that, though not impossible, these scenarios are also unlikely. 

5. Summary and Discussion. 

We have outlined a method that can construct debris trails, using negligible 
computational resources, for any specified satellite mass and orbit and in most parent 
potentials, requiring only the assumption of a timescale over which the satellite is destroyed 
(see §2.2). 

The debris distributions generated by our approach demonstrate a remarkable 
agreement with those calculated from the positions of unbound particles at the end of 
simulations run in a fixed Galactic potential using equivalent satellite and orbit parameters 
(§2.2.3 and Figure ||). This success encourages its use as a tool for understanding the 
evolution of tidal streamers from the accretion of small (m sat /M MW < 1/100) satellites after 
the potential fluctuations due to the initial formation of the Galaxy have died down, and 
for which the limitations mentioned in the following paragraph do not apply. 

The approximations inherent in this prescription (see §2.3) imply that it is not 
well-suited to generating debris distributions if the mass loss rate is nearly constant (e.g. if 
the satellite is on a near circular orbit) or if the satellite's orbit is evolving due to dynamical 
friction with the parent galaxy. Both of these effects could conceivably be taken into 
account with a more sophisticated implementation of our technique. The method is clearly 
inapplicable to problems where the parent's potential is either triaxial or fluctuating as the 
result of a recent merger, such that the debris might be expected to diverge further from 
the satellite's orbit. 

I thank Roeland van der Marel and John Bahcall for invaluable comments on the 
manuscript, Lars Hernquist and Mike Bolte for the collaboration that provided the 
inspiration (and simulations) for this project, and Steve Majewski for many useful 
discussions. The simulations used here were performed at the Pittsburgh Supercomputing 
Center. This work was supported by the Monell Foundation and the Institute for Advanced 
Study SNS Membership Fund. 
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A. Simulations 

The following is a brief sketch of the method and parameter ranges for the simulations, 
which are described fully in |JHB. 



In each simulation we represent the Milky Way by a rigid, three-component potential. 
The disk is described by a Miyamoto-Nagai potential ( Miyamoto fc Nagai 1975*1) , the 



spheroid by a Hernquist (1990) potential and the halo by a logarithmic potential: 

GM disk 

^disk — ■ 



^spher 



R 2 + (a + Vz 2 + b 2 ) 
GM sph 

er 



) 
2 



$haio = v 2 h J Hr 2 + d 2 ). (Al) 

We take M disk = 1.0 x 10 11 , M spher = 3.4 x 10 10 , v halo = 128, a = 6.5, b = 0.26, c = 0.7, and 
d = 12.0, where masses are in M , velocities are in km/s and lengths are in kpc. 

Each satellite is modeled with a collection of 10 4 self-gravitating particles, initially 
represented by a Plummer ( |1911| ) model 

Gm , . „. 

--^T? (A2) 

where m is the mass of the satellite and r is its scale length. During the simulations, the 
mutual interactions of the satellite particles are calculated using a self-consistent field code 
( |Hcrnquist fc (Jstriker 1992 ), with the Milky Way's influence calculated from the equation 



In this paper, we consider the state of the satellite after 10 Gyrs of evolution in twelve 
different simulations, whose parameters are outlined in Table [I]. The first four models 
(1-4) have satellite and orbit parameters chosen at random to explore a range of possible 
outcomes. The remaining eight models all had the same orbit as Model 2, but with satellite 
parameters chosen to change the rate of mass loss and debris dispersal. In Models 5-8 
the satellites have the same central density as Model 2, but different velocity dispersions. 
In Models 9-12 the satellites have the same velocity dispersion as Model 2, but different 
central densities. 
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Fig. 1. — The upper panel plots azimuthal time period as a function of energy for orbits 
with angular momentum J/J c i rc = 0.05 — 1.0 in the spherical Milky Way potential $mw 
(eq. [ |A1|| ). The lower panel plots the ratio of radial to azimuthal time periods Tr/T^ as a 
function of energy. 
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Fig. 2. — The fractional error in the calculation of e (solid line - see eq. |§) and (dotted 
line) if our three-component Milky Way potential (see Appendix A) is approximated as 
purely logarithmic at all radii. 
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Fig. 3.— Snapshot of particle positions in the orbital plane of a simulation of a tidal 
encounter. The dots label particles still bound to the satellite, and the "-" and "+" signs label 
those unbound on orbits with lower and higher energy than the satellite orbit respectively. 
The dashed line shows the satellite's path. The bold circle shows the physical scale calculated 
from equation (|3]). 
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Fig. 4. — Distribution of debris in scaled energy q = AE/e at the end of each simulation 
(histograms). The bold lines show the simple analytic functions that we have used to 
approximate each of these distributions (see eq. ]T2| ). 
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Fig. 5. — Final distribution of debris along orbits for each simulation (histograms). The 
bold lines are the predictions from equation (15). 
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Fig. 6. — Ratio of the density of a remnant following the complete destruction of a satellite 
to the critical density for self-gravity to be important. The different curves are for orbits with 
the same energy and J/ J C j rc = 0.28, 0.36. ..1.0 or R a /R p ~ 0.02 — 1.0. Pericenter corresponds 
to the dip in each curve. 
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Fig. 7. — Mass surface density distribution for debris along the orbits of those satellites with 
known proper motions, assuming each has lost half its mass in the last 10 Gyrs. The dotted 
lines show the current position of each satellite. In the right hand panels, the gap in the 
surface density at this point is the region that is occupied by particles that are still bound 
to the satellite Note that \P is measured from the Galactocentric, not heliocentric position. 
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Fig. 8. — Conversion of Figure [7] to a comparison of the density of horizontal branch stars, 
./Vhb; with the background star counts from the Bahcall-Soneira model for the Galaxy, Nbs- 
The squares in these panels indicate where the orbit is within 20 degrees of the Galactic disk 
and where the signature of the debris is most likely to be swamped by foreground objects. 
The dotted lines show the current position of each satellite. In the right hand panels, the 
gap in the surface density at this point is the region that is occupied by particles that are 
still bound to the satellite 
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Fig. 9. — Evolution of the central surface density of a remnant following the complete 
destruction of a satellite, for orbits with the same energy and eccentricities as in Figure |6] 
Pericenter corresponds to the dip in each curve. 
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Fig. 10.— The time taken for debris to disperse over a full Great Circle (in units of 
the azimuthal time period of the orbit), as a function of distance and for a variety of 
satellite masses (solid lines). The dotted line is the curve T 27r = 10 Gyears. The Galactic 
satellites are plotted either at their current distance or their calculated pericentric distance 
(for those with proper motions). The equivalent time period for a circular orbit at each 
distance is shown along the top axis. For example, this plot implies that it will take 
(T 27r /T^ ~ 19) x ~ 2.3 Gyrs) = 43.7 Gyrs for the tidal streamers from Draco to encircle 
the Galaxy. 
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Fig. 11. — Each curve shows the fractional sky coverage of debris from a satellite of mass 
m sat after 10 Gyrs, as a function of its Galactocentric distance. The dotted lines show the 
behavior if the streamer is not assumed to overlap itself when it extends more than 2-7T along 
the satellite's orbit. 




Fig. 12. — The curves show the number of trails observed along a random tine of sight (or, 
equivalently, the fraction of the sky covered) for total accreted mass M acc as a function of 
the mass of the accreted satellites. 



